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The phase diffusion of the order parameter of trapped Bose-Einstein condensates at temperature 
ksT ^> huj is determined, which gives the fundamental limit of the line-width Av of an atom laser. 
In addition a prediction of the number fluctuations in the condensate and their correlation time r c 
is made and a general relation for Av r c is derived from the fluctuation-dissipation relation. 
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Bose-Einstein condensation in a weakly interacting 
Bose-gas in three dimensions in the thermodynamic limit 
of an infinitely extended system is a second order phase 
transition in which an order parameter, the macroscopic 
wave- function, appears spontaneously with a fixed but 
arbitrary phase, turning the global [/(l)-symmetry re- 
sponsible for particle-number conservation into a broken 
or hidden symmetry [0. The rigidity of the phase of 
the order parameter against local perturbations and the 
absence of any phase diffusion gives rise to the Gold- 
stone modes, which take the form of collisionless (zero)- 
sound or hydrodynamic sound, respectively, depending 
on whether the sound frequency is in the collisionless 
mean-field regime or in the collision-dominated regime 

In finite systems, and thus also in all trapped Bosc- 
gases, sharp phase-transitions are impossible and hidden 
symmetries in a rigorous sense cannot appear Nev- 
ertheless a macroscopic wave-function describing a Bose- 
Einstein condensate still exists, as is now firmly estab- 
lished by the experiments S. However, for the general 
reason mentioned, the phase of the macroscopic wave- 
function cannot be stable but must undergo a diffusion 
process, which restores the U{\) symmetry over suffi- 
ciently long time intervals |EJ. This diffusion process is 
therefore different from the Goldstone modes mentioned 
before, which are oscillations around a fixed value of the 
phase and do not restore the symmetry |J , [|| . 

Recently a first attempt has been made to measure 
the stability of the phase of the macroscopic wave func- 
tion in a trapped BEC. In an experimental set-up of con- 
siderable ingenuity [Q the relative phase of two BEC's 
was measured using a time-domain separated-oscillatory- 
field condensate interferometer. Over the time-interval of 
100 ms scanned in the experiment the relative phase was 
found to be robust. At first sight this experimental re- 
sult may seem surprising since decoherence of entangled 
states of many atoms should be extremely rapid. Then, 
however, one realizes that there is so far no clear theoret- 
ical prediction of the decoherence-time of Bose-Einstein 
condensates in traps against which the aforementioned 



experiment, or extensions of it which will surely follow, 
could be checked. In a number of papers [^) the disper- 
sion of the phase of a trapped Bose-Einstein condensate 
at zero-temperature was considered, which is due to ther- 
modynamic fluctuations Sfi of the chemical potential /i 
in a finite system with fixed particle number. An exten- 
sion of this mechanism to finite temperature has also been 
proposed ||. This effect is not a 'phase diffusion' but cor- 
responds to an effect of inhomogeneous broadening, and 
is even reversible in 'revivals'. The experiments are done 
at 'high' temperature ksT 3> Tiui and even ksT /i, 
where u) is the geometrical mean of the three main trap 
frequencies. One would expect a proper phase-diffusion 
process to occur in such a regime due to the interac- 
tion of the condensate with a thermal bath of collective 
modes and quasiparticles but so far an understanding of 
this process seems lacking (see however This gap in 
our understanding of Bose-Einstein condensates in traps 
at finite temperature is particularly painful, because the 
fundamental limit of the line-width of an atom laser de- 
pends on it: there is not yet a 'Schawlow-Townes' formula 
H for the linewidth of an atom laser, because its deriva- 
tion requires a prior understanding of phase-diffusion in 
Bose-Einstein condensates. 

In the present paper I intend to outline a theory 
of dissipation and thermal fluctuations in a trapped 
Bose-Einstein condensate which is used to determine the 
phase-diffusion constant, and from it, the line- width of a 
trapped Bose-Einstein condensate as a function of tem- 
perature. The result obtained explains the experimen- 
tally observed robustness of the phase. I find it con- 
venient to present first a phcnomcnological framework 
for the theory, in the form of Langevin equations in 
which dissipation appears via phenomenological param- 
eters and the fluctuation-dissipation relation is invoked 
to determine the fluctuations. Then the phenomenolog- 
ical parameters are fixed by drawing on known micro- 
scopic results for the damping of collective modes and 
calculating the new transport coefficient in the Langevin 
equation of the condensate. It is found to result primar- 
ily from scattering of thermally excited collective modes 
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(phonons) off the condensate. This coefficient then deter- 
mines the phase-diffusion constant and the fundamental 
line- width of an atom laser via a Schawlow-Townes-type 
formula. 

The weakly interacting Bose-gas in a trap in standard 
notation is described by the Hamiltonian 



H = I d 3 x 



V 2 +V(x)- (1) 



2m 



The total number of atoms N is fixed, fi is the chemical 
potential. The presence of a Bose-Einstein condensate 
means that many (No 3> 1) particles occupy the normal- 
ized single-particle state ipo (x of lowest energy /i satisfy- 
ing § -(h 2 /2m)V 2 ^ + (V(x) + U No\^o\ 2 )^o = l*Po- 

The number density of the condensate is uq(x) = 
No\ipo(x)\ 2 . Both ipo(x) and No are functions of [i. In the 
following we turn this around and consider \i a function 
of -/V . The presence of the highly occupied condensate 
mode makes the decomposition of the Heisenberg field- 
operator ip(x,t) = (ao(t)ipo( x ) + w(x, i))e~ vt / ri useful, 
where we shall follow the tradition starting with Bogoli- 
ubov [Q and describe the condensate classically, remem- 
bering however, that iVo = |ao| 2 is the particle number 
in the condensate. ip'(x,t) is the field operator for the 
particles outside the condensate. The Hamiltonian then 
splits according to H — H0+H1+H2+H3+H4 where H n 
(n = 0, 1, 2, 3, 4) comprises the terms of H which are of n- 
th order in ip','>p + , respectively. The equation of motion 
for the condensate amplitude ao receives contributions 
from Ho to H3. For a discussion of various microscopic 
approximation-schemes see |To) ] . 

Quasiparticle operators a u ,a+ 

are defined by the standard Bogoliubov transformation 
tp'(x,t) = Y,v {u v (x)a v (t) + v*(x)a+(t)), where u v , v v 
satisfy the usual Bogoliubov-De Gennes equations [fTo) , 
and the a v obey the Heisenberg equations of motion a„ = 
j-[H (u \a u ] with ffM = Tuo v (a+a v - j\v v (x)\ 2 d 3 x). 



Again |u„| 



and uj u are functions of ao 



In order to derive a phenomenological equation of mo- 
tion for the condensate alone we turn to the free energy 
-F(|ao| 2 ). Expanded to second order around its minimum 
it takes the form 



PF(\a \ 2 ) 



(|q | 2 -(K | 2 )) 2 
2<«5A 2 ) 



with (|q | 2 > = N - (N'), (£A 2 ) = (N> ) - (A') 2 , where 

N' = J d 3 xip' (x)ip'(x). The expectation values (N') 
and (SN 2 ) (which turns out to be anomalously large 
~ A 4 / 3 ) have recently been evaluated within the Bogoli- 
ubov theory p^ ] and can therefore here be considered as 
known. The equation of motion of ao near thermal equi- 
librium can now be written with the help of F^ao) 2 ) i n 
the general form JllJ 



ihdo = (1 — £To) 



dF 



Fo(t) 



(2) 



with Gaussian white noise Fo(t) satisfying 

(F (t)) =0, (F*(t)F Q (0)) = 2hk B TT S(t) 

determined so as to ensure the correct equilibrium distri- 
bution p(a ,o%) = Z^ 1 CKp(-F(\ao\ 2 )/k B T) for the 
condensate. Here I ignore the possibility of the occur- 
rence of squeezing in the thermal bath of uncondensed 
particles. Then only a single new phenomenological co- 
efficient To, a dimensionless but surely temperature de- 
pendent number, remains to be determined below. 

As a short digression let us also extend the equations 
of motion of the quasi-particles to include dissipation and 
fluctuation within a phenomenological Markoffian frame- 
work fll5[ . It is convenient to do this by writing quantum 
Langevin equations um 



i%a v = (1 — iT v 



dot 



Fu(t) 



(3) 



with Gaussian Langevin-force operators satisfying 
(F u (t)) = 0, ([F u (t),F+(t')}) = 2h 2 T l/ uj u 8(t~t')5 lJU , and, 
by the fluctuation dissipation theorem, 



(F+(t)F v ,(t')) = 2ti 2 T v u v n v 5(t-t')8 vv , 



(4) 



with the Planck-distribution n v = (e^ hu> " - I)- 1 . The 
phenomenological coefficients r„ have the meaning of one 
half of the inverse Q-factor of mode 1/, r„ = (2Q„) _1 , and 
remain to be determined below. 

But let us now consider how number fluctuations and 
the phase-diffusion of the condensate in equilibrium fol- 
lows from (||). Its deterministic part describes the re- 
laxation of the condensate to the minimum of the free 
energy F of the condensate at (Ao) = (|ao| 2 ). The parti- 
cle number fluctuations 5No = \ao\ 2 ~~ (|ao| 2 ) m equilib- 
rium, after linearizing (||) in 5 No, are found to have the 
correlation function 



(8N (t)SNo(t')) = (5AT 2 ) e -l t - t 'l/^ 



(5) 



with the correlation time t c = (h(5N§) /2To(N }k B T) 
which could in principle be measured by taking the 
Fourier transformation in time of time-resolved in situ 
phase-contrast images of the condensate. On a time- 
scale very much larger than the correlation time t c the 
phase fo of the condensate in equilibrium, i.e. the phase 
of ao — "v/Aoe^ , satisfies the Langevin equation of a 
Wiener process with diffusion constant 



D v = k B T(To + T : 1 )/(h(No}), 



(6) 



i.e (((f (t) — ^o(O)) ) = D v t. The expectation 
value (ao(t)) then decays exponentially according to 
(a (t)) = J(No)e-^ vt with the linewidth Av given by 
the Schawlow-Townes-type formula 

Av = k B T(T Q + T Q 1 )/(2h(N )) > k B T/(hN Q ) . (7) 
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The general relation between Av and r c 



The relevant matrix elements are 



Av 



h(SN 2 ) 



2h(N )^2(N }k B TT c 



2(N )k B TT c 



(8) 



follows from our phenomenological theory, which is inde- 
pendent of the yet unknown coefficient T and holds for 
general value of k B T (outside the critical region). Eq. (||) 
is a general consequence of the fluctuation dissipation re- 
lation for the condensate in the absence of squeezing in 
the bath of uncondensed atoms. 

Let us now determine the numbers To, r„ from mi- 
croscopic considerations, starting with Y v for v 7^ 0. We 
shall here confine our attention to the damping of the low- 
lying collective modes in the collisionless regime, even 
though our phenomenological framework may still be 
used in the collision-dominated regime. Furthermore, we 
confine ourselves to the experimentally relevant regime 
k B T ^> fi. The damping of the low lying collective 
modes in this regime is due to Landau damping, de- 
scribed by a part of H3, and was calculated for spa- 
tially homogeneous condensates first by Szepfalusy and 
Kondor [0. Their result written for our coeficient T v 
yields Y v = A u (fc B T/ J u)(n a 3 ) 1 ^ 2 - The numerical coef- 
ficient A v turns out to be independent of v in the spa- 
tially homogeneous system, and its value there is 
A v = Sit 3 / 2 /A. However, a similar expression for T u was 
even shown to hold for the collective excitations in traps 
|jl9f where no is the condensate density no(0) in the cen- 
ter of the trap and A v is a numerical coefficient which 
depends on the trap geometry and the mode function for 
mode v. For the dipole modes A v must vanish by the 
Kohn-theorem po| . 

The coefficient Tq is unknown so far and needs to 
be calculated from scratch. Here we shall determine 
the temperature dependence of Tq for large quasi- 
homogeneous condensates for which the local density ap- 
proximation is applicable, collecting the trap-dependence 
of Tq in a dimensionless prefactor which we leave un- 
determined. Collisions of quasiparticles with the con- 
densate changing the particle number |ao| 2 in the con- 
densate by A|ao| 2 = ±1 are described by H 3 = 
J d 3 xUo(aotpo(' l P + ) 2 ip' + h.c). In such processes the en- 
ergy changes only by a tiny amount AF = dF/d\a \ 2 
(of the order of k B T j \J (SN'q)). However, processes de- 
scribed by H3 involve three quasiparticles besides one 
condensate particle, and they can therefore take up, 
from a slightly perturbed condensate, an arbitrarily small 
amount AF of its free energy. 

By the golden rule the rate 7 = d(\ao\ 2 )/dt is given by 



2tt 



+ AF/h) 



~ \M™M -o>v-t»*- AF/R)) (9) 
[n v %(n K + I) - (n„ + l)(n p + l)n K ] 



M£j,n = 2U a J d 3 xip v v (u%u^ + ^v^v^) + (v <-> \j) 
M^ K = 2U a a J d 3 xtp ut(v*V K + |tt*u«) + {v <-> /i) 

AfW describes a Landau-scattering process in which one 
atom is scattered out of the condensate by the absorp- 
tion of the two quasiparticles v, fi out of and the emission 
of the new quasiparticle k into the thermal bath. Like- 
wise describes Beliaev scattering where an incoming 
thermal quasiparticle k is absorbed, an atom is kicked out 
of the condensate, and two quasiparticles v, [i are emitted 
into the thermal bath. The factor 1/2 in (0) accounts for 
the indistinguishability of pairs v 1 fi and /i, v. 

In the phonon part of the excitation spectrum we have 

—1/2 

u\ ~ — v\ ~ u> x . Furthermore, in that low- energy 
region the statistical factor [...] in (j|) is well approxi- 
mated by (k B T) 2 (uj K — ujjy — Lo^/f^LOjyLo^uj^, where the 
frequency difference in the nominator becomes ±AF/h in 
the product with the 8- functions, which express energy- 
conservation. Anywhere else the small energy AF is neg- 
ligible. The frequency factors in the denominator, to- 
gether with similar factors in the denominator coming 
from the matrix elements, make the phonon contribu- 
tion to the sums in Q the dominant one, at least in 
large condensates, and we shall therefore concentrate on 
this contribution in the following. This frequency range 
has a natural upper cut-off at fj,/h, where the collective 
phonons go over smoothly into particle-like excitations. 
In finite condensates also a natural lower cut-off exists 
at about the trap frequency Q, where the phonon wave- 
length becomes comparable with the size of the conden- 
sate. 

We shall evaluate the sums in (||) in the usual Thomas- 
Fermi and local density approximation plj ] for u\ , v\ , u>\ , 
integrating over the wave vectors of the phonons in an 
interval implied by the cut-offs. We obtain in this way 



7= {2n)3h8 (AF|oo| ) 



dwudio,, 



W v U),j,(u v + Up) 



Clearly the contribution near the lower cut-off at Q dom- 
inates in the double frequency integral. After its evalu- 
ation, and using Uq = Airh 2 a/m and do — ijfi/muj to 
introduce the s-wave scattering length a and the zero- 
point amplitude do , we obtain to leading order in fiui / /i 
d{\a \ 2 )/dt = 7 = -2(T /h)((dF/d\a \ 2 )\a \ 2 ) with 



To = A, 



18 In 2 ,ak B T,2 
7T dohui 



(10) 



Here A$ depends on the trap geometry and is of order 
1 in an isotropic trap. For the experimentally realized 
condensates r << 1 is implied by (|To|) . 

In order to give a practically useful estimate of the 
phase-diffusion rate let us compare the line-width Av of 
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the condensate with the directly measurable inverse life- 
time r" 1 of a collective mode k with frequency of order 
Q. We find for T << 1 



[4] D. S. Hall, M. R. M atthews, C. E. W ieman, and 



Avt k 



Uj/UJ K 



(»/k B Tf 



A A K 2881n2 L (A^ )K(0)a 3 ) 3 / 2 



(11) 



For n/k B T = HT 1 , (N ) = 10 6 ,n (0)a 3 = 3.HT 6 the 
factor [...] on the right hand side of ( pi] ) is about 2. This 
may explain the experimentally observed Q robustness of 
the phase. Another way to put our result is to note that 
up to a numerical prefactor the relation fuDTo /// ~ T% = 
(2Q K )~ 2 holds. Hence, according to (fj]) for r << 1 



Av = (A 2 JA 36ln2)(k B Tu>Q 2 Jn(N )) 



(12) 



where Q K is known experimentally. 

In conclusion we have provided a theory of thermally 
generated phase-diffusion in Bose-Einstein condensates 
and determined a temperature dependent Schawlow- 
Townes type formula for the linewidth Av of the con- 
densate. This is also the minimum line-width of an atom 
laser based on the Bose-Einstein condensate. We have 
also calculated the correlation time r c of the particle 
number fluctuations in the condensate which should be 
measurable. Last not least we have derived a relation 
(^J) between this correlation time r c and the line-width 
Av which follows from the fluctuation-dissipation rela- 
tion independently of any microscopic detail, but subject 
to the assumption of negligible squeezing in the thermal 
bath seen by the condensate. Removing this restriction 
will be the subject of a more detailed paper. 
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